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Abstract 

We obtained new upper limits on the diffuse gamma rays from the inner Galactic 
(IG) and outer Galactic (OG) planes in 3-10 TeV region, using the Tibet air shower 
data and new Monte Carlo simulation results. A difference of the effective area of 
the air-shower array for observing gamma rays and cosmic rays was carefully taken 
into account in this analysis, resulting in that the flux upper limits of the diffuse 
TeV gamma rays were reduced by factors of 4.0~3.7 for 3~10 TeV than those in our 
previous results (Amenomori et al., 2002, ApJ, 580, 887). This new result suggests 
that the inverse power index of the energy spectrum of source electrons responsible 
for generating diffuse TeV gamma rays through inverse Compton effect should be 
steeper than 2.2 and 2.1 for IG and OG planes, respectively, with 99%C.L.. 
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1 Introduction 



The EGRET observation (|Hunter et al.l . 119971 ) of diffuse gamma rays up to 



10 GeV from the inner Galactic (IG) and outer Galactic (OG) planes shows a 
sharp ridge along the Galactic plane. The flux in the lower energy region ( E < 
1 GeV) is consistent with the conventional calculation (e.g. 'DermeJ,'l986f), al- 



though the EGRET flux is higher than the COS B data (Mayer-Hasselwander et al, 
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1980l ) by a factor of about 3 in several GeV. In order to explain the EGRET 
flux in terms of cosmic-ray hadronic interactions with the interst ellar raatter 
(ISM), a flatter source proton spectrum with power index of 2.45 (jMoril . 119971 ) 
or 2.25 (I Webber . 19991 ) is required ne ar the Galactic center with th e thickness 
|6| < 10° or < 5°. On the other hand, Porter and Protherod (1997) suggested 
that cosmic-ray electro ns, in such high energies, may play a significant role 



Pohl 



sucn mgn energies, may play a signmcant role 
and EsDositol 11998') showed that if the source 



in gamma-ray creation, 
electron has a power index of - 2.0, the EGRET excess above 1 GeV can 
be well explained by the inverse Compton (IC). Also, Berezinskv et al. ( 19931 ) 
calculated the diffuse gamma-ray flux in the wide energy range of TeV-PeV 
based on the cosmic-ray hadronic interactions with ISM. 

In the energy region above 1 TeV, almost all observations have given upper 
limits on the flux of diffuse gamma rays from the Galactic plane, except for 



the recent MILAGRO result giving a 4.5 a excess (jPlevsher et al.L 1200,^^ at 
1 TeV for 40° < I < 100°, |6| < 5°, though their inner Galactic plane in- 
cludes the Cygnus region where some complex features are observed in wave- 
length bands of radio, in frared and GeV gamma rays. In the previous paper 
^Ameno mori et al. . 20021 ) we presented the flux upper limits of diffuse gamma 
rays at 3 TeV and 10 TeV from IG and OG planes using the Tibet-II and 
Tibet-Ill data, assuming that the effective area of the array is the same for 
air sho wers initiated by gamma rays and cosmic rays with the same primary 
energy (jAmenomori et al.L l200lh . Since then, the Tibet air shower array has 
been gradually enlarged and improved by adding more detectors, covering an 
area of 37000 m^ in which 789 detector are placed at the grid of 7.5 m spac- 
ing. We have also done detailed Monte Carlo (MC) simulations to examine 
the performance of the ne w air shower array for o bserving cosmic rays and 
gamma rays above 1 TeV ( Amenomori et a'lll2003| ). 



In this paper, we present the results on the flux of diffuse gamma rays ob 
tained with the Tibet a ir shower array, based on the n ew MC simulation 
(lAmenomori et aL , 2003[ ) and compare the previous results ( Amenomori et al 
20021) . 



2 Monte Carlo simulation and data analysis 



A MC simulation for performance of the Tibet array (Pig. 1) is a key subject 
in this analysis. However, we briefly describe its outline because of a page 
restrictions. In this simulation, primary cosmic rays were thrown isotropically 
into the atmosphere with an appropriate energy spectrum and composition. 
In practice, the heavy dominant (HD) model, given in Pig. 2, was employed, 
which can well des cribe the experimental data in the wide energy rage from 
1 TeV to 10^ TeV (|Amenomori et al.l . bood I2OOIL l2003l ) . The primary gamma 
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rays were assumed to arrive at the Earth with an inverse power spectrum with 
indices of /?=2.2, 2.4, 2.6, and 2.8 both from the I G plane (20° < I < 55°J &| < 
2°) and OG plane (140° < / < 225°, |6| < 2°) (|Amenomori et all 120021 ). as 
shown later in Fig. 6. 

Here, it may be worthwhile to note that at a high altitude like Yangbajing 
(606 g/cm^) the size of air showers induced by gamma rays is larger than that 
induced by protons and other nuclei in the multi-TeV region. The shower size 
ratio of gamma-ray origin versus proton origin takes a value around 2 ~ 3 
for the same primary energy. Hence, the detection efficiency of gamma rays is 
better than that of protons at Yangbajing. 

Simulated events were generated with the same array configuration as the 
experiment shown in Fig. 1 and they were analyzed under the same condi- 
tions and procedures as the experimental data. It was confirmed that the 
mode energies, ^'mode, of proton- induced showers observed by the Tibet-Ill 
and Tibet-II are 3 TeV and 10 TeV, respectively, where the Tibet-II means a 
part of the Tibet-Ill overlapping with a detector spacing of 15 m. The Tibet 
II and HI data were analyzed in 4° width binning of the IG and OG planes 
in the equatorial coordinates as shown in Fig. 3, because the respective 120 
curved belts are quite the same both in the zenith angle distribution and in 
the average primary energy of arriving air showers from the IG and OG planes, 
respectively. 



: FT Detector 
: FT Detector W/D-PMT 
Density Detector 



Fig. 1. Layout of the Tibet III array in 2000. 

The effective areas of the Tibet HI array for gamma rays, protons and cosmic 
rays from the direction of IG plane are compared at the given energy bins in 
Fig. 4. It can be seen that the effective area for gamma rays is about 7 times 
larger than for the galactic cosmic rays in the multi-TeV region, as already 



reported elsewhere (jAmenomori et al.l . |2 
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Fig. 2. Primary cosmic-ray spectrum in the HD model. 
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Fig. 3. Binning for the IG and OG planes in the analysis. 
The energy spectra of simulated events coming from the sky region of IG 
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Fig. 4. Effective areas of the Tibet III for IG plane. 

plane, for example, are shown in the differential form in Fig. 5, where the 
spectral index of 2.6 is assumed for both gamma rays and cosmic rays. We can 
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see that the mode energy of triggered gamma rays is 1.5~2.0 times smaller 
than the galactic cosmic rays for the same trigger condition. Hence, if this 
difference is taken into account, the average advantage factor of the effective 
area for gamma rays is estimated to reduce to 4.0 for i?mode —3 TeV and 3.7 
for 10 TeV in average of the gamma-ray spectral indices of (3 = 2.2 ~ 2.8. 
In the present analysis, the median value (3=2.5 is chosen because of its weak 
dependence on the spectral index. 
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Fig. 5. Primary spectra and distribution of triggered events for IG plane. 



Array of data taken with 


Tibet III 


Tibet II 


-^mode 


3 TeV 


10 TeV 


Inner or Outer Galactic plane 


I G 


G 


I G 


G 


Significance of excess (cr)^ 


+2.52 


+0.25 


+1.71 


-0.63 


Flux ratio of gamma rays versus cosmic rays^^ 
/^(at 1 a)/IcR = 1/VB (xlO-4) 


1.95 


1.16 


2.43 


1.45 


99%C.L. upper limit in the original (/3=2.4)^ 
E?^dN{> E^)/dE^ (xlO-3cm-2s-isr-iMeV) 


9.6 


3.3 


4.0 


1.3 


Effective area (7/CR) ratio of the Tibet array 

5cff(7)/5'eff(CR) 


4.0 


4.0 


3.7 


3.7 


99%C.L. revised upper limit with /9=2.5 
E^dN{> E^)/dE^ (xlO-3cm-2s-isr-iMeV) 


2.6 


0.88 


1.2 


0.37 



Table 1 

Effective area ratios and gamma-ray flux upper li mits. 

The rows with ^ are referred to the original paper ( Amenomori et al. . 20021 ) . 
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It is also found that the effective area ratio gives almost non-difference on the 
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Fig. 6. Gamma-ray spectra drawn from 10 TeV. 
sky region between IG and OG planes. 



3 Results and discussions 



In this experiment, no significant signal of diffuse TeV gamma ray s at the 



Earth was detected from the Galactic planes. As in the previous paper (jAmenomori et al 



2m± 



the significance (a) of the excess events of TeV gamma rays from the 



IG or OG plane was calculated using a simple formula of {E — B)/y B, where 
E is the number of events on the on-plane and B the number of background 
events estimated from the neighboring bins around the on-plane. In Table 1, 
those evaluated values by the previous analysis are tabulated together with 
the simulated effective area ratio of gamma rays versus cosmic rays and the 
revised upper limits with a small change of the gamma-ray spectral index (3 
from 2.4 to 2.5. In this table, the 3 TeV data were obtained by the Tibet III 
array (inner area of 22,050 m^), and the 10 TeV data by the Tibet II array 
(28,350 m2). 

The flux upper limits, thus obtained, may be compared with the EGRET 
data. As an example, we show the differential spectrum of diffuse gamma rays 
from the IG plane in Fig. 6. The straight lines denote the interpolation with 
the assumed spectral indices of 2.2, 2.4, 2.6 and 2.8 at the Earth, drawn from 
the upper limits at 10 TeV. From this flgure, we may say that the index 2.5 is 
the most probable value in the GeV- TeV region, if the considerable portion of 
the EGRET bump in th e 1~10 GeV region is possibly attributed to n° decay 
( Aharonian et al.l . l2000h . 
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In Figs. 7 and 8, we plot the new flux upper limits at -E'mode=3 TeV (Ti- 
bet III array) and 10 Te V (Tibet II array) together with the previous ones 
f Amenomori et al. . 2002f l. New data give more stringent upper limits than the 
previous ones and this difference can be mostly attributed to a detailed calcu- 
lation of the effective area of the array for the observation of ga mma rays and 
cosmic rays. Also, shown in these flgures are the EGRE T data ( Hunter et al.1 . 



1997t ) of \b\ < 2°, the upper li mits by Whipple fW) (|LeBohec et al.L 



200C 



with 99.9%C.L., HEGRA (H) (|Aharonian et all 1200 if) with 99%C.L., where 
the 2nd and the 3rd experiments scanned a sma ll sky region arou i id I = 40°, 
and also shown are HEGRA- AIROBIC C (Ha) (|Aharonian et all l2002h and 
CASA-MIA (C-A) (|Borione et al.L Il998l ) both with 90%C.L. for co r nparis on. 
For the only positive result from MILAGRO (M) (|Flevsher et al.L l200,5l ). it 
should be noted that the observation of signiflcant excess is located at the 
Cygnus region rather than in the IG plane according to their signiflcance map. 



The calculations for the inverse Compto n (IC) were done by Porter and Protherod 
ffl997^ rPP2.0 and PP2.4 in flgures), and lTatevama and Nishimurai ((2003i)(TN2.0 
and TN2.4 in flgures), where 2.0 and 2.4 are source electron spectral indices. 
The theoretic al curve for gamma ray s of hadronic origin through 7r° decay is 
also given bv lBerezinskv et all (|l993l )rBGHS). 
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Fig. 7. Diffuse gamma rays from the IG plane of 20° < I < 55°, \ b\ < 2°. 



In conclusion, our new upper limits with 99%C.L. can give a constraint that 
the spectral indices of source electrons for IC should be steeper than 2.2 in the 
IG plane and also 2.1 in the OG plane in comparison with the IC theoretical 
curves, when the power index of 2.5 is adequate for the arriving diffuse gamma 
rays in the multi-TeV region. 
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Outer Galactic plane 




ENERGY (GeV) 

Fig. 8. Diffuse gamma rays from the OG plane of 140° < / < 225°, \b\ < 2°. 
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